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Abstract 

We have explored spectral line profiles due to spiral patterns in accretion discs around black holes. A 
parametrization was employed for the shape and emissivity of spiral waves, which can be produced by 
non-axisymmetric perturbations affecting the disc density and ionization structure. The effects of the light- 
travel time, energy shift, and gravitational focusing near to a rotating black hole were taken into account. 
A high-resolution ray-tracing code was used to follow the time variations of the synthetic line profile. A 
variety of expected spectral features were examined and the scheme applied to a broad iron line observed 
in MCG-6-30-15. 

Key words: accretion, accretion disks - black hole physics - galaxies: nuclei - galaxies: individual 
(MCG-6-30-15) - X-rays: galaxies 
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Introduction 

Time-dependent phenomena are of extreme richness, 
variety, and observational relevance for accreting black- 
hole systems (e.g., Kato et al. 1998; Krohk 1999). 
Clearly, astrophysically realistic processes must result 
from the interplay of different factors which affect light- 
curves and produce spectral features. Here we concen- 
trate on the X-ray spectral band. 

Rapid temporal changes of the flux and of individ- 
ual spectral features are observed in many active galac- 
tic nuclei (AGN) and in Galactic black-hole candidates 
(GBHCs). Variability patterns fit well with the widely 
adopted scenario of these objects as massive black holes 
surrounded by accretion flows, in the form of a disc or a 
torus (e.g., Ulrich et al. 1997). The X-ray variability of 
AGN and of GBHCs is often attributed to instabilities 
in their accretion flows causing changes in the inner and 
outer boundaries of the radiating disc region or relatively 
persistent non-axisymmetric inhomogeneities in the disc 
itself. This paper deals with variability properties which 
can be due to spiral structures forming in the inner part 
of the disc at distances < lOOrg from the centre. We 
consider rotating (Kerr) black holes together with spiral 
waves of density perturbation in the disc matter. The 

1 Low-resolution figures have been used in order to eonform to 
the standards of astro-ph e-print Hbrary. 



waves are supposed to modify the emissivity pattern of 
the gaseous flow. 

A power-law continuum is observed in X-ray spectra, 
which is thought to provide primary irradiation of the 
flow. Numerical computations of local reflection spectra 
from X-ray irradiated material suggest that the emissiv- 
ity depends strongly on the nature and location of the 
primary source (e.g., Martocchia, Matt 1996; Martocchia 
2000, and references therein), and also on the ionization 
structure and density profile of the flow (Done, Nayak- 
shin 2001). Only in the case of density larger than a 
critical value and small ionization parameter, the matter 
behaves like a pure optically thick medium from which a 
typical cold iron line is emitted (Nayakshin et al. 2000; 
Dumont et al. 2000). The extent to which ionization in 
the disk influences emission and absorption features is 
subject of much debate. 

As far as the power-law continuum is concerned, a 
plausible scheme suggests that this spectral component is 
produced in a hot Comptonizing medium (corona) above 
the disc. Another possibility is that magnetic reconnec- 
tion causes temporary, strong flashes of radiation (flares) 
near the inflowing matter (Galeev et al. 1979; Coppi 
1992; Haardt et al. 1994; Nayakshin 1998). Several mech- 
anisms have been proposed for flares, but the available 
spectroscopy does not allow us to decide among them; 
see Malzac and Jourdain (2000) for recent discussion and 
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further references. The flares are supposed to irradiate 

the gas in their neighbourhood, contributing to short- 
term variabihty (Mineshige et al. 1994; Yonehara et al. 
1997). Furthermore, the coronal activity may be local- 
ized in asymmetric and rapidly evolving areas. In both 
cases (flares vs. inhomogeneous coronal activity) varia- 
tions in the primary flux are expected, which also de- 
termine the variation in the additional spectral features 
due to Compton reflection (a continuum "hump") and 
fluorescence (lines). 

In this work we explored the expected spectral line pro- 
files in the region where relativistic light-travel time, en- 
ergy shift and lensing effects influence substantially the 
observed radiation fluxes, temporal variability, and the 
lines centroid energy. We assumed that a perturbation 
(of density and ionization structure) develops on length- 
scales of ten or a few tens of gravitational radii, giving 
rise to emissivity contrast against the background of the 
axisymmetric flow. The spirals become active after being 
illuminated from the primary source. Here, one has to 
distinguish between two cases: (i) the primary illumina- 
tion independent from the origin of the spiral perturba- 
tion of the disc; (ii) the two phenomena (the disc pertur- 
bation and the flare) being connected to each other. The 
second case is particularly relevant for binary systems 
with Roche-lobe overflow; in such a situation, the outer 
end of the spiral pattern coincides with the place where 
the overflowing material hits the disc. This is the place 
where a flare or a spot develops. Therefore, we consider 
the emissivity variations with simultaneous contributions 
from the spiral pattern and from the flare. 

We computed the predicted, variable line profiles with 
the aim of assessing the importance of the mentioned ef- 
fects, especially the degree of the line flux and energy 
fluctuations that could be assigned to the existence of 
such non-axial structures. We also mention our compu- 
tational approach, which can be adapted to different in- 
trinsic emissivitics and the source geometry, and, for the 
sake of illustration, we determine free parameters of the 
phenomenological model for the variable iron Ka emis- 
sion of the Scyfcrt galaxy MCG-6-30-15. 

Rapid variations have been reported in the X-ray con- 
tinuum and in the Ka iron line of MCG— 6-30-15 (6- 
7 keV), namely, Iwasawa et al. (1996, 1999) described the 
spectral changes which occurred during ASCA observa- 
tions in 1994 and 1997. It is important to examine how 
the changes of selected spectral features — and the flu- 
orescent iron feature in particular — reflect the changes 
in primary illumination, because this helps to constrain 
the geometry of the emitting regions, for instance by re- 
verberation techniques (Reynolds et al. 1999; for a recent 
review, see Fabian et al. 2000). 

We start the next section by giving some further hints 
on the physical models which justify our interest in con- 
sidering the spiral-type patterns. 



2. The Model 

In AGN, variability time-scales go down to less than 
one hour, which is almost comparable with the light- 
crossing time (<c) at distances of the order of the grav- 
itational radius (rg) of a supermassive black hole: tc ~ 
10^M7 s, where M7 is the central mass in units of lO^M©. 
Other time-scales relevant to black-hole accretion discs 
(orbital, thermal, sound-crossing, and viscous) are typ- 
ically longer than tc- Characteristic time-scales of the 
observed variability also depend on the inclination of the 
system with respect to an observer, and can be scaled 
with the mass of the central object. Moreover, time inter- 
vals become longer very near to the hole if gravitational 
delays predicted by general relativity are taken into con- 
sideration. 

2.1. Astrophysical Motivations 

The irradiation of the disc from a varying source of 
primary X-rays docs not lead to instantaneous response, 
but generates reflection or fluorescence. A delay occurs 
depending on the geometrical and physical state of the 
gaseous material. For cixamplc, Blandford and McKcc 
(1982) and Taylor (1996) have discussed the intrinsic de- 
lays of the response in the case of clouds geometry, but 
here such effects arc not considered for simplicity. In- 
stead, it is the location of the primary source — often 
identifled in one or more flaring regions — together with 
reflection patterns in the disc which determine the ob- 
served signal. A single variability event is therefore made 
of a primary flare and a complex response due to Comp- 
ton reflection (or fluorescence, in the case of the iron line) 
by the disc matter at various distances from the hole. 

Several people have considered the reprocessing of pri- 
mary X-rays recently (for references, sec Diimont et al. 
2000; Karas et al. 2000; Nayakshin et al. 2000). Repro- 
cessed radiation reaches the observer from different re- 
gions of the system. Individual rays experience unequal 
time lags of 6t for purely geometrical reasons, and other 
relativistic effects occur in the very vicinity of the hole 
and are relevant for the source variability, namely, lens- 
ing and Doppler boosting (Ranch, Blandford 1994). In 
the case of long-lived patterns (t » tc) in the disc, there 
is, indeed, a substantial contribution to the variability 
caused by the orbital motion (Abramowicz et al. 1991; 
Mangalam, Wiita 1993). If the disc is strictly equato- 
rial, two types of such patterns have been examined in 
more detail: spots (which could be identified with vor- 
tices in gaseous discs; cf. Abramowicz et al. 1992; Adams, 
Watkins 1995; Karas 1997), and spiral waves. In princi- 
ple it is possible to distinguish different origins of the 
varying emission, namely, short-lived flares which decay 
on the dynamical time-scale versus long-duration features 
(Lawrence, Papadakis 1993; Kawaguchi et al. 2000). To 
that purpose, the geometry of the system needs to be con- 
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strained sufficiently. In addition to intrinsic variations, 
obscuration is involved along the path of light rays to- 
wards the observer. The light propagation is affected by 
the presence of refractive interstellar medium on the line 
of sight (Rickett et al. 1984) which contributes, together 
with internal and random obscuration, to rapid and high- 
amplitude flickering seen by the distant observer (Abras- 
sart, Czerny 2000). These non- intrinsic effects are be- 
yond the scope of this paper. 

The spiral waves represent large-scale structures (size 
comparable with the radius), which can be induced by 
non-axisymmetric perturbations (Sanders, Huntley 1976; 
Wada 1994; Lee, Goodman 2000) or due to magnetic in- 
stabilities (Tagger et al. 1990). Also, a pattern resem- 
bling a single-armed spiral is produced from an extended 
spot after its decay due to shear motion in the differen- 
tially rotating disc. Even though such a pattern is only 
a transient feature, it may last sufhciently long to pro- 
duce observable effects. Thus, the spirals are expected to 
arise naturally in accretion discs and to persist there for 
times longer than the orbital period of the flow particles. 
Sanders, Huntley (1976) demonstrated that the forma- 
tion of a two-arm spiral wave is a natural response of dif- 
ferentially rotating gaseous discs to an oval gravitational 
perturbation, much like in the case of a self-gravitating 
stellar system into which the gas is assumed to be em- 
bedded. Recently, Kuo and Yuan (1999) examined the 
excitation of spiral waves in Lindblad resonances by com- 
paring the asymptotic theory of bar-driven waves with 
the corresponding results of numerical computations. In 
a different way. Tagger and Pellat (1999) showed that 
spirals can be formed due to a certain kind of magnetic 
instability; such waves remain localized close to the inner 
radius of the disc. In general, spiral patterns can develop 
an odd or even number of arms depending on the type of 
the perturbing forces that act on the disc. 

Spiral patterns have been treated in the theory of stel- 
lar systems, and widely discussed in the context of non- 
axisymmetric perturbations which can modify density 
profiles and emissivities of gas discs (Goldreich, Tremaine 
1979; Sawada et al. 1986; Rozyczka, Spruit 1993). The 
status of SPH computations was summarized recently 
by Lanzafame, Maravigna and Belvedere (2000) who ex- 
amined the role of stellar-mass ratio on the formation of 
spiral structures in binaries. The spirals meet numerous 
applications also with respect to cataclysmic variables 
(Steeghs, Stehle 1999). Furthermore, the modulation of 
X-rays by spiral waves was also proposed with regard to 
polars (Murray et al. 1999). 

Most of the attention with regard to spiral- waves spec- 
tral features has been motivated by studies of cataclysmic 
variables and other stellar-mass binary systems. It was 
noted that the variation in the density profile and ioniza- 
tion of accretion flows, predicted by numerical and semi- 
analytical methods, is followed by temperature modu- 
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Fig. 1. Geometry of the model: an artistic expression of 
spiral patterns in an accretion disc. In the present 
work, the intrinsic emissivity is described by a small 
number of free parameters which also determine the 
variability of the spectral features. We consider the 
emissivity to be modulated by the spiral pattern 
with an associated decaying flare. 



lation, and therefore by a change in the gas (thermal) 
emissivity within the spiral patterns. A similar effect 
is also expected for the X-ray irradiated accretion flows 
in AGN due to an increase in the matter density along 
the patterns (if optical thickness t < 1) or variations in 
the ionization structure (the emissivity is not sensitive to 
density changes if r ^ 1). A phcnomcnological descrip- 
tion appears to be adequate for the current scenario of 
X-rays reprocessing in AGN cores, and is convenient for 
the numerical computation of line profiles. 

As far as AGNs are concerned, spiral waves have been 
explored as possible sources of spectral line variations in 
the context of broad emission-line regions (Chen, Halpern 
1990; Chakrabarti, Wiita 1993, 1994; Eracleous 1998) 
and in maser sources (Maoz, McKee 1998). In these stud- 
ies, because the radiation originates at distances > 10'^7-g 
from the centre, general-relativity effects are unimpor- 
tant and variability time-scales are of the order of several 
months to years. On the other hand, Sanbuichi, Fukue 
and Kojima (1994) considered spirals extending close to 
a Schwarzschild black hole where the effects of general 
relativity play a role, and showed examples of relativisti- 
cally distorted spectra. 

2. 2. A Toy Model of Spiral- Wave Emissivity 

We parametrize the Kerr spacetime in terms of the 
black hole specific angular momentum (0 < a < 1). In 
the disc plane, we consider both the spiral and a decay- 
ing spot locahzed at the spiral's outer end (radial loca- 
tion and azimuth of the spot are thus determined by the 
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Fig. 2. Maximum time lag Stm shown between the rays 
from a ring of given radius in the Kerr metric (solid 
curves). Four cases are plotted for radii greater than 
rms and different combinations with the black-hole 
rotation parameter, a. The observer inclination an- 
gle 6»o is: (a) a = I, Bo = 80°; (b) a = 1, do = 20°; 
(c) a = 0, 00 = 80°; (d) a = 0, = 20°. Ra- 
dius is given in units of rg, while the time lag is 
expressed with respect to Keplerian orbital period, 
tit. The Euclidean estimate of <5tm is also shown 
(dotted lines). 



spiral's shape and its orientation with respect to the ob- 
server). Such a spot can represent a temporary ofF-axis 
flare (its maximum amplitude is a free parameter). 
See figure 1 for a general scheme of the model. 

The complexities of primary X-ray reprocessing are 
overpassed by directly parametrizing the emissivity. Our 
approach to the emissivity parametrization is close to the 
computations of Sanbuichi et al. (1994) and Bao and Wi- 
ita (1999), who considered a non-rotating black hole. 

The spiral pattern is characterized by the emissivity 
distribution in the equatorial plane, 



Vo(?-c), 



(1) 



where the frequency-dependent term, jo, is taken as a 
narrow Gaussian profile, while the power-law index, 7, 
describes the overall radial decrease in the emissivity. 
The spiral rotates with respect to the observer, and 
its shape is logarithmic and determined by the function 
(Pe{r) = if + Q!log(r/ro), where arctana is the pitch an- 
gle, and ro is the outer radius where the wave is being 
excited.^ This parametrization approximates the spiral 

2 The locally emitted intensity (denoted by subscript "e" ) is re- 
lated to the observed intensity at spatial infinity (subscript "o") 
by the formula Je(i'e; t — 5t)/u^ = jo(i'o; t)/i>^, where the ratio 
g = Uo/ye defines the redshift factor and St = 5t{rc,ipe',6o,o,) 
is time delay due to the light-travel time along the ray. See 
Fanton et al. (1997) and Martocchia, Karas and Matt (2000) 
for further details on ray-tracing computations. 



pattern evolving on the background of a Keplerian disc. 
The azimuthal form of the spiral depends also on the 
type of perturbation, so what we use here is one of the 
possible modes, which appears to be simple enough for 
illustration purposes. 

As already discussed, the model takes inspiration from 
spirals in close binaries where a spot arises following 
the interaction between the stream (after the Roche-lobe 
overflow) and the disc (e.g., Wolf et al. 1998). The radia- 
tion flux in the spot decays with an e-folding time of g > 1 

orbital periods, = 27r(ro^^ -I- a). Again, the frequency- 
dependent part, Jo, of the spot intrinsic profile is taken in 
the Gaussian form and is assumed to be sufficiently nar- 
row (its smearing is due mainly to radial component of 
the source orbital motion), while the amplitude, A{r, ip), 
decreases exponentially with the distance from the flare 
center located at rg. 



(2) 



Here, A{r,ip) — yl„iexp(— Ar) with, typically, the spot 
size Ar = |i — tqI ^ 1 and ro ~ 20. Hence, the resulting 
emissivity is a sum of contributions (0) and (|^): jc = 

As) , (f) 

Jc + Jc ■ 

It is the exact form of the emissivity distribution along 
the disc surface which determines the observed flux in the 
line wings and, in particular, which of the two wings — 
the bluish or the reddish one, appears more pronounced 
in the predicted spectrum. Thereby the lines are more 
complex than one would expect on the basis of intu- 
ition with simple ring- type sources (double-horn profiles). 
Here, we assume that a spiral with two arms extends 
from r = vq down to r^s, being terminated at the inner- 
most stable orbit. Single-armed spirals show comparable 
spectral characteristics, but their variations occur with a 
double period and a somewhat higher amplitude, similar 
to the case of an elongated spot. 

2.3. Light-Travel Time Near a Black Hole 

Now we concentrate our attention on the consequences 
of finite light-travel time upon line variations. 

One can write an order-of-magnitude (Euclidean) es- 
timate of the maximum light-travel time delay, Stm, of 
rays originating in an equatorial ring (radius rc centered 

on the black hole): 5tai{rB,9o) ~ 7r~^rc sin 6*0 (6*0 is 
the observer's inclination). Here, St^ characterizes the 
mutual delay of the signals emitted at different parts of 
the ring, and is defined as the span in arrival times of 
the rays originating at different (fc-^ Usually, r^ is sup- 

3 Light rays were approximated as straight lines for the pur- 
pose of estimating Stjn here. For consistency with subsequent 
paragraphs, the standard notation will be adopted for Kerr 
spacetime in Boyer-Lindquist coordinates (Misner et al. 1973). 
Geometrized units c = G = 1 are used and all lengths are made 
dimensionless by expressing them in units of the typical mass 
of the central black hole, M ^ M7 = IO'^Mq . 
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Fig. 3. Isocontours of constant light-travel time St and of redshift function g are shown in equatorial plane of the Kerr black 
hole. The contours are shown for a non-rotating black hole (a = 0, left panels), and for a maximally rotating hole (a = 1, 
right panels). Observer inclination is, from top to bottom, 9o = 20°, 9o = 50°, and 6o = 80°. For a = 0, three circles 
are plotted with radii of the horizon, of the circular photon orbit, and of the marginally stable orbit. The circles coincide 
with each other if a = 1 (corotating orbits are assumed). 
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posed to be greater than a marginally stable orbit, r^s'- 
fms = 3rg for a non-rotating black hole with a = 0, and 
'"ms for a maximally rotating black hole, a — *■ 1 

(rg = 1 + \/l — a^). Light-travel time can be ignored 
provided St^ <C 1. On the other hand, the lags become 
increasingly important when the radiation source extends 
very close to the black-hole horizon where the resulting 
situation is further complicated by the lensing effect. 

The geometrical time lag can be roughly characterized 
by 5tm, which indicates whether the Euclidean formula 
estimates the correct value of the light-travel time with 
an acceptable precision, or whether relativistic effects 
modify the actual delays significantly. Figure 2 gives Stm 
for a source located at > r^ns- The solid curves repre- 
sent the values of the delay computed in Kerr spacetime, 
while the dotted lines (slope —1/2) show the approxima- 
tion. Relativistic corrections become increasingly impor- 
tant for Te/rg < 3, in which case St^ increases sharply. 
On the other hand, the difference between the exact value 
of 6tm and its Euclidean approximation is less than 10 % 
for a source location > 5rg. 

Figure 3 shows the absolute value of light-travel time 
difference St between a reference ray (arbitrarily chosen 
but fixed) and a ray coming from a given point of emis- 
sion {re,(fe} in the equatorial plane {9 = 90°). In this 
figure, St was calculated in Kerr metric (no approxima- 
tion was used here). Given a and 9o, the values of 6t 
determine where in the disc the corrections on the time 
lag play an important role. Reference values of 5t can be 
transformed to physical units, as measured by a distant 
observer, by the relation St [s] ~ lO^My St. The equa- 
torial plane is portrayed from above (along the common 
rotation axis of the hole and of the disc), and a distant 
observer is located at the top (Cartesian y ^ oo; co- 
ordinates X and y are introduced in the disc plane by 
x"^ + y^ = rl). Contours of constant time lag are almost 
horizontal far from the hole. A distortion of the con- 
tours from straight lines occurs near to the hole due to a 
lensing effect and gravitational delay. (Indeed, the con- 
tour lines can even become split into two disjoint parts.) 
Next, the redshift function, g, is also plotted by another 
set of contour lines in figure 3. Relativistic effects on the 
time lag and redshift are more important for a maximally 
rotating black hole (right panels) because of the smaller 

2.4- Parameters of Spiral- Wave Line Profiles 

Spiral-wave structures produce a rich variety of spec- 
tra and light curves, much broader than we can observe 
in the case of localized spots. Apart from the black hole 
angular momentum, a, and the orientation angle of the 
source with respect to observer, 9o, one can identify sev- 
eral parameters related to the geometry and orbital mo- 
tion of such spirals which determine the spectral variabil- 



ity: 

• The non-axisymmetric intrinsic emissivity of the source 
is described by equation (|l]) with free parameters (a, f3, 
7, ro) determining the azimuthal and radial dependence of 
emissivity. 

• Doppler shifts between the emitted and detected radia- 
tion increase the centroid energy and intensity of the rays 
originating at the approaching side of the source, and vice 
versa for the receding side (factor g"^). This effect tends to 
enlarge the width of the spectral features and to enhance 
the observed flux at its high-energy tail. 

• Gravitational redshift affects mainly photons coming 
from the inner parts of the source. This effect decreases 
the centroid energy of the line. 

• Light-travel time determines, along with the shape of 
the spiral, which photons are received at the same instant 
of time. We remark that contours of constant time delay 
are deformed near the hole. Hence, variability amplitude 
can be much accentuated if the spiral coincides with one of 
the contours at certain orbital phase. This effect influences 
the flux from some patterns under suitable orientation. 

• Gravitational lensing enhances the flux from the far side 
of the source (the part of the source which is at the upper 
conjunction with the hole). This tends to increase the ob- 
served flux around the centroid energy. 

All parameters of the intrinsic emissivity and of the 
black hole spacetime enter in computations of the indi- 
vidual ray trajectories which form the resulting radiation 
flux at the detector. We find that the dependence of ob- 
served profiles on 7 is rather weak for a typical choice of 
7 ^ 2. One should be aware of other implicit pa- 
rameters which were kept fixed in the present model: the 
inner radius of the spiral pattern (here r^s), the (small) 
width of the intrinsic Gaussian local line profile, the num- 
ber of spiral patterns (here n = 2) , and duration of the 
flare (> tk)- In all of these computations, for the sake of 
simplicity the radiation was emitted isotropically in the 
rest frame of the gaseous medium. 

3. Results 

A typical line profile of a double-armed spiral is shown 
in figure 4 where the observed spectral feature is plotted 
as a function of the observed energy and orbital phase. 
The variation in the spectrum is shown over the whole 
revolution of the pattern by 27r in azimuth (as the nor- 
malized phase ranges from to 1), while energy axis is 
normalized to the local emission energy, i/g. Contours of 
the dynamical spectrum are shown by projecting the sur- 
face plot onto the plane of energy vs. phase. Here, one 
can clearly identify that the contribution of the spiral 
alone is repeated twice in each of the subplots (because 
we considered double arms), but the total signal does not 
bear this property because of the decaying flare contribu- 
tion. Next, figures 5-6 show these contours for different 
model parameters. The contours are plotted with uni- 
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Fig. 4. Example of the dynamical spectrum of a line from a spiral wave near a black hole (a = 0). Left panel: The predicted flux 
is background subtracted and normalized to maximum. The energy-integrated lightcurve and time-integrated spectrum 
are also shown. The parameters for this plot were a = 6, /3 = 18, 7 = 0, ro = 13, Am = 1. The inclination angle was kept 
fixed at 50°. Right panel: Projection is shown onto normalized energy vs. phase diagram with flux indicated by different 
levels of shading. 



form spacings between the maximum flux and zero flux 
of the signal from which the underlying continuum is sub- 
tracted. 

3.1. The Predicted Form of Observed Profiles 

The parameter space has been explored systematically 
in the following manner. 

• Black-hole angular momentum: a = (columns denoted 
by letters a, c, e, g), and a = 1 (b, d, f, h). The two cases 
correspond to a non-rotating and to maximally rotating 
hole, respectively. 

• Observer inclination: 60 — 20° (the rows denoted by 
roman numbers i-iii), 60 = 50° (iv-vi), and 9^ = 80° (vii- 
ix). The disc plane is 60 = 90°. 

• The outer edge of the spiral pattern: vq = 7 (i, iv, vii), 
ro = 13 (ii, V, viii), and ro = 20 (iii, vi, ix). Units of GM/c? 
are used. 

• The form of spiral pattern, as defined by its pitch angle 
arctana and contrast f) in eq. (|l|): a = 3, /3 = 8 (a, b); 
a = V3, /3 = 8 (c, d); a = 6, /? = 8 (e, f); a = 6, /3 = 18 
(g, h). Large values of (3 correspond to well-defined spirals 
(with small width and large contrast), while large values of 
a define tightly wound spirals (corresponding to high Mach 
number in the disc medium). 

• A contribution of the flare to the line in terms of its 
duration parameter, q — 2 (units of the orbital period), 
and the initial flux. Am. The average ratio of the flare-to- 
spiral intrinsic emissivity in the centre of the line was taken 
as Am = 1:2* for a weak flare (flgure 5) and Am = 1:1 
for a strong one (figure 6). 

In each of the sub-figures the spectra were computed 
with sufficient resolution in energy and phase (10^ x 10^ 



grid) by the ray-tracing code (Martocchia 2000). 

One can recognize several characteristic features of the 
dynamical spectra. The flare is easily distinguished in 
figure 6 where it was assumed to be more intense with 
respect to the spiral pattern than in figure 5. However, 
even in the latter situation the Doppler boosted radiation 
of the flare dominates over the extended spiral when the 
disc is viewed at large inclination, and especially if the 
whole source is near to the horizon (ro ~ 7; see the row 
vii). Next, one can recognize the effect of gravitational 
redshift if the source is seen almost along the rotation 
axis (i-iii), whereas this effect turns out to be less visible 
in very broad spectral features corresponding to large 
inclinations (vii-ix). 

3.2. The Case of MCG- 6-30-15 

While the two plates in figures 5-6 has given us a gen- 
eral feeling for the infiuence of different parameters, a 
suitable (x^ minimizing) routine is used in practice for 
fitting the computed profiles to data. 

Figure 7 shows the best fits to three iron-line profiles 
in MCG— 6-30-15 (we use exactly the same data as de- 
fined by Iwasawa et al. 1996). These correspond to the 
spectral states of the source during its 1994 observation 
with the ASCA satellite. Here, we do not intend to ex- 
actly mimic what is going on in the source, because this 
appears to be premature with the available data. It is, 
however, possible to constrain the distance from the cen- 
tre at which the line could originate, depending on the 
angular momentum of the hole and inclination of the ob- 
server. 
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Fig. 5. Dynamical spectra of the two-arm pattern with a slowly decaying flare. Different cases correspond to the form of the 
spiral, rotation parameter of the black hole, and observer's view angle (see the text for details). Here, the spiral dominates 
over the flare. 
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Fig. 6. Same as in figure 5, but now the flare is assumed to be more pronounced. This represents an intermediate situation 
between a pure spiral and the orbiting spot model. 
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7. Left: Time variation of the best-fit profile dur- 
ing the entire revolution of the spiral. Right: The 
computed profile is plotted after averaging to three 
different phases. The feature observed in MCG— 6- 
30-15 (Iwasawa et al. 1996) is overlayed with an 
appropriate phase shift. 



The curves shown in the right panels are those to which 
the minimizing routine converges in the parameter 
space (as described in subsection The reduced 

is (from top to bottom): xV^-o.f. = 72/53, 40/55, 60/40 
for the three phases, respectively. Here, a = 0, 6o = 20°, 
ro = 13, a = 6, /3 = 18, = 1/16. We thus find that 
the data are well fitted by the double arm at moderate 
inclination angles. In such a case, the dependence of the 
observed profile on a is negligible because a substantial 
part of the line flux originates sufhciently above r^g- 

4. Discussion 

In the present work we did not intend to ameliorate the 
fits of MCG— 6-30-15 data further. It is, however, worth 
noticing that even the adopted, rather simple scheme pro- 
duces acceptable fits and the required, reddened line cen- 
troid also with small values of a. One can thus address 
the question if those redshifted and variable profiles of 
the iron line could be achieved without the necessity of 
invoking a large mass (> 1O''M0) of the central object 
and/or its fast rotation, which have recently been put 
in certain doubt by several people. It is thus legitimate 
to consider all different viable options because various 
effects contribute to the line formation, and the inter- 
pretation of observed profiles is not straightforward.^ A 

4 Weaver and Yaqoob (1998) and Reynolds and Begelman (1997) 
proposed alternative schemes which do not involve fast rotation 
of the black hole. Nowak and Chiang (2000) discussed the 



self-consistent interpretation of the profiles cannot be 
achieved without regard to the underlying continuum 
component where it is crucial to distinguish between line 
and continuum photons. 

Our discussion is relevant for X-ray variability of some 
spectral states observed in AGN and GBHCs; it is quite 
obvious, however, that only a part of the complex vari- 
ability behaviour could be attributed to motion of large 
scale structures, such as spiral waves. The model is com- 
plementary to the discussion of instant off-axis flares in- 
ducing the observed line variations, which enable one to 
restrict the geometry of the source by means of the rever- 
beration technique. Proper evaluation of the light-travel 
time from different parts of the source is important in 
both schemes. 

5. Conclusions 

We computed variable line proflles due to non- 
axisymmetric spiral-type patterns near a Kerr black hole. 
We adopted several simplifying assumptions and con- 
strained some of the parameters to fixed values. Yet, we 
could show that a big variety of line profiles are repro- 
duced and the data of MCG— 6-30-15 which we used for 
illustration purposes can be modeled without the neces- 
sity of a — > 1, although such a possibility is not rejected. 
Because the existence of spiral features is expected on 
the basis of theoretical models, our discussion seems to 
be better substantiated than data fitting by ad hoc func- 
tions. Also, the adopted approach uses computed pro- 
files, rather than being restricted to integral quantities 
(centroid energy, line width, etc.); it can thus provide ro- 
bust results once better-resolved data become available. 

Our present scheme can be considered as comple- 
mentary to the picture of Reynolds et al. (1999) and 
Ruszkowski (2000), who examined off-axis flares above 
the disc surface as the origin of line variations, and to 
Vaughan and Edclson (2001), who addressed the problem 
of the link between line and continuum variations. Natu- 
rally, the phenomenological description of the source will 
have to be supplemented by a time-dependent physical 
model of X-ray reprocessing. 

Although we applied our computations only to the 
case of MCG— 6-30-15, there is a more general inter- 
est in the problem of spectral features from inner parts 
of gaseous discs endowed with spirals. Compact binaries 
with matter overflow represent another class of objects 
where our computations could be useful. Patterns like 
those discussed in this paper can be introduced in or- 
der to explain asymmetric lines and variability observed 

possible caveats in interpretations of the present time-resolved 
spectra, and they concluded that the central mass is about 
10® Mq, much less than previously considered. A rather small 
mass, however, requires efficient emission, which also seems to 
pose a problem, so the issue has not been settled yet. 
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in UV/optical spectral bands from the broad- line region 

(Dumont, CoUin-Soufrin 1990; Corbin 1997), and peri- 
odicities reported in the optical light curve of the blazar 
OJ287 (Takalo 1994; Villata et al. 1998). In these cases, 
however, the radiation is emitted at distances exceeding 
10^ gravitational radii, so that the general relativity cal- 
culations are unnecessary and the computations can be 
performed in a much faster way. 

Even though the results cannot be conclusive with the 
present-day resolution, there is evident potentiality of 
this approach with planned missions, like Constellation- 
X and XEUS. The model of spirals illustrates the op- 
portunity of tracing the surface structures on accretion 
discs. 

We acknowledge support from the grants GACR 
205/00/1685 and 202/99/0261. VK thanks for hospi- 
tality of the International School for Advanced Studies 
(Trieste) and Observatoire de Paris-Meudon. 
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